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Abstract 

Polymer  electrolyte  fuel  cells  (PEFC)  in  the  regime  of  limiting  current  is  considered.  Maximum  current  limitations  come  due  to  oxygen 
deplete  on  the  cathode.  They  appear  strongly  dependent  on  the  geometry  of  the  gas  diffuser  and  supply  channels.  Effect  of  liquid  water  in 
the  electrode  and  air  pressure  in  gas  channels  is  studied  to  identify  the  best  geometry  of  the  gas  supply  net  in  the  terms  of  the  maximum 
cell  current.  Comparison  between  model  predictions  and  experimental  data  is  carried  out. 

©  2003  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Polymer  electrolyte  fuel  cells  (PEFC)  represent  increas¬ 
ingly  popular  alternative  sources  of  energy.  However,  fuel 
cells  with  high  power  densities  must  be  developed  for  this 
technology  to  become  a  viable  power  source  for  transporta¬ 
tion  applications.  Three  main  factors,  limiting  the  current 
and  leading  to  voltage  losses  in  proton  conducting  membrane 
fuel  cells  are  the  following:  activation  resistance,  attributed 
to  the  catalyst  layer  in  contact  with  electrode  and  accessible 
by  reacting  gases  [1-4],  ohmic  voltage  losses,  attributed  to 
the  electronic,  ionic  and  contact  resistance  of  the  cell  [5-9] 
and  the  mass  transport  resistance,  when  the  reactant  gases 
deplete  on  the  reaction  interface  as  their  transport  to  the  re¬ 
action  sites  fails  to  keep  up  with  the  reaction  rate.  During  the 
operation,  reactant  gases  flow  down  the  channels,  diffusing 
into  the  electrodes  through  the  diffusion  layers  to  the  cat¬ 
alyst  interfaces.  At  high  current  densities,  the  reaction  rate 
becomes  mass-transfer  limited  by  the  transport  of  the  reac¬ 
tant  gases  from  the  channel  to  the  catalyst  surface.  The  last 
phenomenon,  present  in  particular  on  the  cathode,  still  has 
not  been  fully  addressed  by  researchers.  The  only  common 
method  to  overcome  such  transport  limitations  and  to  in¬ 
crease  the  power  density  has  been  to  use  pressured  reactants. 

So  far,  numerous  theoretical  models  of  PEFC  have  been 
developed  (see,  for  example,  [10-17]).  However,  they  were 
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usually  restricted  to  simple  one-dimensional  cases.  Recently, 
some  researchers  introduced  full  three  dimensional  models, 
accounting  even  for  the  water  transport  in  porous  matter 
[18-29].  Recent  experimental  studies  demonstrated  signifi¬ 
cant  influence  of  final  pressure  drops  and  geometry  of  sup¬ 
ply  channels  on  the  effectiveness  of  the  cell  [28-31].  For 
example,  interdigitated  and  serpentine  configurations  of  sup¬ 
ply  channels  essentially  make  use  of  a  new  pressure  driven 
mechanism  of  mass  transport  that  acts  additionally  to  com¬ 
mon  molecular  diffusion. 

In  the  present  study  we  restrict  our  attention  to  the  oxygen 
electrode,  which  has  been  the  subject  of  numerous  studies 
and  which  is  the  limiting  electrode,  compared  to  the  rela¬ 
tively  efficient  hydrogen  electrode  of  PEFC.  At  the  cathode, 
the  reaction  rate  is  limited  by  the  transport  of  oxygen  from 
the  channels  to  the  catalyst  sides  and  of  the  water  vapour  in 
the  opposite  direction.  In  addition  to  the  gas  transport  lim¬ 
itations  the  cathode  faces  another  problem,  called  electrode 
flooding.  The  water  generation  rate  at  the  cathode  at  high 
current  densities  often  exceeds  its  removal  rate  from  the  in¬ 
ner  porous  layers  in  the  form  of  water  vapour,  resulting  in 
condensation. 

The  gas  forced  flow  from  the  gas  channel  into  the  dif¬ 
fuser,  add  convection  mechanism  of  mass  transport  to  pure 
diffusion.  Since  convection  is  much  faster  than  diffusion,  the 
reaction  rates  at  the  catalyst  can  be  significantly  improved. 
Furthermore,  the  shear  force  of  this  gas  flow  field  helps  re¬ 
move  most  of  the  liquid  water  that  is  entrapped  into  the  inner 
layers  of  the  electrode,  reducing  the  flooding  problem.  This 
effect  is  present  not  only  in  so  called  interdigitated  channels 
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with  closed  ends,  where  gas  is  forced  to  flow  through  the 
electrode  to  exit,  but  also  in  some  other  geometric  configu¬ 
rations,  such  as  serpentine  channels.  However,  the  air  meets 
high  flow  resistance  in  interdigitated  flow  and  it  results  in 
low  airflow  rates,  low  outlet  oxygen  concentration.  On  the 
contrary,  serpentine  design  makes  use  of  both  advantages 
of  the  conventional  and  interdigitated  flow  fields.  When  ser¬ 
pentine  channels  in  the  electrode  are  set  under  pressure,  the 
flow  rate  of  air  along  the  channels  becomes  sufficiently  high, 
providing  enough  oxygen  supply  for  the  device,  whereas 
oxygen  is  transported  to  the  inner  layers  of  the  electrode  by 
both  diffusion  and  convection  that  is  induced  by  the  pressure 
drop  between  the  neighbouring  channels. 

We  first  report  the  experimental  part,  then  we  describe  the 
two-dimensional  electrode  model  mathematically  and  apply 
the  model  for  the  comparative  analysis  of  the  maximum  cur¬ 
rent  limitations  of  the  PEFC,  running  it  along  the  3D  di¬ 
mension  for  specific  configurations  of  gas  supply  channels 
of  serpentine  design.  In  particular,  the  nets  of  2,  3,  4,  5,  7 
and  10  parallel  serpentine  channels  on  a  square  plate  will 
be  considered.  The  aim  of  the  studies  is  the  identification 
of  most  efficient  in  the  terms  of  the  limiting  current,  chan¬ 
nel  configuration.  Then,  the  model  results  are  analysed  and 
compared  with  experimental  results. 

2.  Experimental 

2.7.  Materials 

Commercially  available  20  wt.%  Pt/C  catalyst  powders 
on  carbon  black  (Vulcan  XC72)  were  purchased  from  E-Tek 
Inc.  Three-layer  (substrate/diffusive  layer/catalyst  layer) 
gas  diffusion  anode  and  cathode  (106  cm2)  were  prepared 


using  a  spray  technique  described  in  detail  in  previous 
works  [32-34].  The  substrate  was  carbon  paper  (Toray 
TGPH090).  The  weight  composition  of  the  diffusion  layer 
was  85  wt.%  of  carbon  and  15  wt.%  of  PTFE,  with  carbon 
loading  of  1.93  mg  cm-2.  The  catalyst  layer  was  prepared 
by  mixing  appropriate  amounts  of  carbon- supported  cata¬ 
lyst  (24.5  wt.%)  and  5  wt.%  Nation  ionomer  solution  from 
Aldrich  (14.9  wt.%)  and  glycerol  (60.5  wt.%).  The  plat¬ 
inum  loading  in  all  anodes  and  cathodes  was  kept  constant 
at  0.34  mg  cm-2  on  the  MEA.  Nation  115  membrane  (Du 
Pont)  was  used  after  purification  treatment  in  5%  (w/v) 
H2O2  solution  at  80  °C  for  1  h,  followed  by  a  second  treat¬ 
ment  in  1  M  H2SO4.  The  ME  As  were  formed  by  hot  press¬ 
ing  the  electrodes  (106  cm2)  onto  the  membrane  at  130  °C 
for  1-5  min  and  50-100  kg  cm-2. 

Membrane  electrode  gasket  assembly  (MEGA)  technol¬ 
ogy  [32]  was  used  and  a  well-defined  shape  compatible  with 
the  cell  hardware  was  achieved  (Fig.  l).This  technology  al¬ 
lows  to  disassemble  the  cell  and  to  replace  the  graphite 
plates,  maintaining  the  same  MEA  [32].  In  this  work,  alu¬ 
minium  anticorrodal  100  (Alusuisse)  end  plates  (185  mm  x 
185  mm  x  11  mm)  were  used.  At  the  lateral  borders  of  the 
assembly  and  end  plates,  manifolds  holes  for  gas  feeding  are 
visible  (Fig.  1).  Graphite  flow  field  plates  were  assembled 
with  a  typical  parallel  channel  configuration  (3  or  10  chan¬ 
nels  on  the  cathode  side  and  10  channels  on  the  anode  side) 
using  commercially  available  BMA5  graphite  produced  by 
SGL  Carbon  Group  (Germany). 

2.2.  Electrochemical  characterisation 

Full-cell  electrochemical  tests  were  carried  out  in  the  10 
or  three  channels  systems  using  a  106  cm2  single  cell  in¬ 
corporated  in  a  Globe  Tech  Inc.  mod.  890  test  station.  Two 


Fig.  1.  Schematic  of  a  single  PEM  fuel  cell  showing  the  location  of  the  components:  (A)  aluminium  or  SS316L  end  plates,  (B)  copper  current 
collector/gaskets,  (C)  graphite  flow  field  plates,  (D)  MEGA. 
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Fig.  2.  Cell  voltage  and  power  density  vs.  current  density  for  three 
(squares)  and  10  (circles)  cathode  channels  configuration,  at  Tce\\  = 
75  °C  under  H2  14  see  min-1  A-1  and  O2  52.5  see  min-1  A-1  flux  at 
1.5/1. 5  absbar. 


aluminium  end  plates  and  two  graphite  current  collectors 
composed  the  single  cell.  The  MEGA  system  and  the  other 
components  of  single  fuel  cell  are  schematised  in  Fig.  1. 

Cell  voltage  versus  current  density  measurements  were 
galvanostatically  performed  with  a  programmable  power 
supply  interfaced  with  a  computer  for  data  acquisition.  The 
results  are  given  in  Fig.  2. 

All  measurements  were  carried  at  75  °C  in  H2/O2  flux 
at  1.5/1. 5  bar  absolute  pressure.  The  temperature  humidifier 
was  kept  at  85  °C  for  the  anode  and  50  °C  for  the  cathode. 


3.  Description  of  the  model 

A  porous  electrode  is  a  very  complicated  structure,  which 
is  difficult  to  characterise  on  a  microscopic  level.  One  ap¬ 
proach  to  this  problem  is  a  porous  electrode  theory  [35].  It 
employs  the  macroscopic  model  that  accounts  for  the  essen¬ 
tial  features  of  a  porous  electrode  without  going  deep  into 
exact  details.  The  development  of  the  equations  for  describ¬ 
ing  the  membrane  oxygen  transport  in  the  gas  diffuser  in  this 
model  is  based  on  such  macrohomogenious  description  with 
uniform  morphological  properties  such  as  porosity,  tortuos¬ 
ity,  and  permeability.  In  actual  fuel  cells  and  especially  in 
multi-cell  stacks,  thermal  gradients  are  likely  to  apply  due  to 
the  heats  of  reaction  on  the  electrodes  and  the  product  water 
is  likely  to  be  a  mixture  of  liquid  and  vapour.  However,  with 
the  available  apparatus  it  is  hardly  possible  to  characterise 
these  gradients  and  state  of  water  precisely.  Moreover,  the 
temperature  gradient  in  the  cathode  effects  density  and  dif- 
fusivity  coefficients,  but  the  density  diffusivity  product  that 
characterises  the  oxygen  transport  varies  less  than  4%  for 
practical  temperatures  in  PEFC  if  the  entering  the  channel 
air  is  fully  saturated  [28]. 

3.1.  Assumptions 

In  order  to  analyse  the  diffusion  effects  on  cathode  be¬ 
haviour  and  address  current  limitations,  we  assumed  several 
simplified  approximations  for  system,  which  are  commonly 
made  in  the  current  literature.  These  are  as  follows: 


•  All  the  porous  zones  in  the  fuel  cell  domain  are  assumed 
to  be  isotropic  with  uniform  macrocharacteristics  such  as 
porosity,  tortuosity,  and  permeability — liquid  water  mod¬ 
ifies  the  effective  porosity  for  the  oxygen  gas  transport  in 
the  cathode  diffusion  layer. 

•  We  assume  isothermal  properties  for  the  system  and  ne¬ 
glect  the  effect  of  gravity. 

•  The  gas  mixture  is  considered  a  perfect  gas. 

•  In  the  pores  of  the  gas  diffuser,  we  assume  that  the  water 
vapour  is  equilibrated  with  liquid  water  in  the  form  of 
small  water  drops.  Though  different  approaches  have  been 
adopted  by  different  researches  [36,37],  in  the  present 
work  we  assume  that  liquid  water  does  not  constitute  a 
continuous  barrier  to  oxygen  transport  [38]. 

•  Water  leaves  the  channels  and  enters  the  electrode  in  the 
form  of  saturated  water  vapour  only. 

•  We  also  assume  that  the  water  in  the  gas  diffuser  takes 
off  only  a  negligible  amount  of  dissolved  oxygen  and 
nitrogen,  compared  to  the  large  gas  phase. 

•  The  active  catalyst  layer  is  assumed  to  exist  only  as  an 
ultra  thin  plane,  so  that  the  transport  of  the  reactants  within 
this  layer  and  the  gas  diffusion  in  the  membrane  can  be 
neglected. 

•  Usually,  gas  pressures  are  approximately  equal  to  each 
other  at  the  anode  and  the  cathode  of  PEFC.  Moreover,  the 
membrane  is  practically  impermeable  for  the  gas  phases. 
Then  the  assumption  that  the  pressure  varies  in  the  dif¬ 
fusers  plane  only — p(x)  =  const  is  justified. 

•  The  processes  in  fuel  cell  are  steady  state  and  the  condi¬ 
tions  are  stationary. 

•  We  also  assume  that  there  are  no  sinks  and  sources  in 
the  body  of  the  diffuser  and  the  electrochemical  reaction 
takes  place  on  the  surface  of  the  catalyst  layer  only. 

Other  assumptions  will  be  discussed  as  the  model  equa¬ 
tions  are  presented  and  the  model  is  applied  to  the  electrode 
of  the  polymer  electrolyte  fuel  cell  with  serpentine  gas  chan¬ 
nel  net.  How  the  liquid  water  is  transferred  in  the  cathode 
is  not  included  in  the  model. 

3.2.  Governing  equations 

First,  consider  two  parallel  gas  channels,  which  are  distant 
at  d  under  a  diffuser  of  thickness  h.  The  schematic  drawing 
of  the  electrode  of  the  polymer  electrolyte  fuel  cell  is  given 
in  Fig.  3. 

Then,  the  fluxes  of  gases  in  the  body  of  the  diffuser  are 
governed  by  well-known  equations  of  diffusion,  continuity 
equations  and  D’Arcy  flow  equation  [39]: 


h  =  -pDiVcA  +  gpiu,  V//  =  0,  u  = - V  p, 

p 

Y  C,  =  1 ,  c;  =  — ,  p  =  pRT,  (1) 

P 


where  7Z  ,  are  the  fluxes  of  i  substances;  c\  are  their  concentra¬ 
tions;  pi  are  their  densities;  p  is  total  density;  p  is  pressure; 


98 


K.  Zhukovsky,  A.  Pozio  /  Journal  of  Power  Sources  130  (2004)  95-105 


▲ 


outlet  gas 
supply  < 
channel, 
pressure  p2 


cA(0,d)=c0 

|Bx(°-y)=0 

lAx(°-y)=0 


Solid 

non  porous 
wall 


cA(0,0)=c 


.0 


inlet  gas 
supply 
channel,  - 
pressure  p1 


02,N2  gas 

mixture 

▲  ▲ 


P(d)  =  P; 


catalyst 

surface 


diffusion 


porous 

gas 

diffuser 

layer 

1 


"*■  Oxygen 
J  flux 


D'Arcy 

flow 


0 


cA(h,y)=o, 
Oxygen  is 
fully 

consumed 


02,N2  gas 
mixture 


p(0)=p1  h 


Fig.  3.  Schematic  drawing  of  an  electrode  of  a  polymer  electrolyte  fuel 
cell. 


Di  are  diffusion  coefficients  of  gases  into  the  mixture;  u  is  the 
bulk  speed  of  the  gas  mixture,  R  is  gas  constant;  T  is  temper¬ 
ature;  k  is  permeability  coefficient;  /c  is  viscosity  coefficient. 

In  order  to  address  oxygen  transport  limitations,  oxy¬ 
gen  of  initial  concentration  co  is  considered  to  be  fully 
converted  on  the  catalyst  surface  of  PEFC.  Boundary  con¬ 
ditions  in  such  case  of  limiting  current  read  as  follows 
(see  Fig.  3): 


In psy  =  -5.80022  x  103/T  +  1.39150 

-4.86402  x  10_2(7’  + 4.176  x  1(T5)7’2 
-  1.44521  x  10~J(T3  +  6.54597)ln  T.  (4) 


The  binary  diffusion  coefficients  are  taken  from  Slattery 
and  Bird  estimation  [41],  which  at  low  temperatures  shows 
that  they  are  dependent  on  critical  temperatures  Tcj,  Tcj 
pressures  pcy,  pcj  and  masses  M/,  Mj  of  the  respective 
components  i  and  j ,  pressure  of  gas  p ,  and  temperature  T,  as 
represented  by  equation 


pDi-j  =  a 


x 


b 

(Pc,iPcj)1/3(Tc4Tcjp12 


( 1  1 Y/2 

\M,  +  Mj) 


where  a  =  3.640  x  10  4,  b  =  2.334  for  water-oxygen 
and  water-nitrogen  and  a  =  3.64  x  10-4,  b  =  1.823  for 
oxygen-nitrogen  diffusion  coefficients. 

With  oxygen  concentration  varying  from  zero  in  catalyst 
layer  to  the  maximum  possible  value  in  saturated  humid  air, 
the  change  in  density-diffusivity  product  pDo2  is  negligible 
(less  than  4%)  under  normal  fuel  cell  operation  conditions 
[28].  Then  at  low  pressures  and  constant  temperature  the 
following  combinations  of  variables  are  constants 


pD  o2 
RT 


=  const, 


=  const. 


Taking  into  account  the  expressions  (3)-(6),  relating  pres¬ 
sure,  temperature,  diffusion  coefficient,  permeability  and 
other  parameters  of  the  gas  mixture  and  porous  diffuser,  the 
above  set  of  Eq.  (1)  with  boundary  conditions  (2)  possesses 
simple  analytical  solution  for  the  concentration  cq2  (. x ,  y) 


P  =  p(y)p(0)  =  pu  p(d)  =  p2,  I,,x(x  =  0,  }>)  =  0, 
c02(h,  y)  =  0,  cq2(0,  0)  =  c0,  cO2(0,  d)  =  c\.  (2) 


Diffusion  coefficient  Dt  is  a  function  of  temperature,  pres¬ 
sure  and  gas  composition.  The  effective  oxygen  diffusion 
coefficient  in  nitrogen-water  vapour  mixture  can  be  calcu¬ 
lated  as  a  function  of  binary  diffusion  coefficients 


1  —  (co2Mo2/M) 

(<W Dq2-w)  +  (cn2/Do2-N2)  ’ 


Ow 


P 


YJCiMi  =  M, 


where,  cw  is  the  concentration  of  water  vapour,  co2 — the 
concentration  of  oxygen,  Do2-w — diffusion  coefficients  of 
oxygen  in  water  vapour,  Do2-n2 — diffusion  coefficients  of 
oxygen  in  nitrogen,  M — the  molecular  mass  of  the  mixture, 
Mi — the  molecular  mass  of  i  substance;  other  notations  are 
the  same  as  in  (1).  Water  saturation  pressure  is  calculated 
from  the  following  semiempirical  approximation  [40] : 


cq2  (x,y)  =  cos  Fx  ^ 


exp 


J  '  '  Z, 


Co  -  Cl  exp  {— ((G  +  H)/2)d } 

1  —  exp  {— Gd) 

—co  exp  {—Gd}  +  ci  exp  {— ((G  +  H)/2)d } 

1  —  exp  {—Gd} 


(7a) 


G  =  J  A  F2  +  H2,  F=—,  H=--3Z — ^1. 

2  h  a  2d 

(7b) 

Here  p\,  P2  are  the  pressures  in  the  two  points  of  the  two 
channels  in  question,  d  is  the  distance  between  those  points, 
h  is  the  thickness  of  the  diffusive  layer,  co — initial  concen¬ 
tration  of  oxygen  in  “in”  channel;  ci — final  concentration  of 
oxygen  in  “out”  channel;  a  and  /3  are  given  by  formulae  (6). 
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3.3.  Solution  for  the  oxygen  flux  in  the  diffuser  between 
two  parallel  channels 

The  expressions  (7),  (7a)  and  (7b)  yield  the  oxygen  flux 
through  the  membrane  Jo2  as  follows: 

ed 

Jo2=  lAx(h,y)dy 

Jo 

71  (  El  H-G  E2  H  +  G  \ 
=  a  —  - exp - y  H - exp - y  . 

h  \H  —  G  F  2  7  H  +  G  F  2  y) 

(8) 

The  above  developed  simple  analytical  model  predicts 
that  for  the  case  of  two  parallel  channels  at  y  =  0  and 
y  =  d,  the  oxygen  flux  per  unit  of  channel  length  through 
the  catalyst  surface  Jo2  is  at  its  maximum,  determined  by  a 
transcendental  function  of  the  ratio  q  =  d/h  of  the  distance 
between  the  two  channels  to  the  diffuser  thickness,  as  well 
as  of  the  physical  parameters  ratio  alp  and  the  difference 
P\  ~  Pi'  Diffuser  with  permeability  k  ~  10-7  to  10-8  cm2 
and  p\  =  105  Pa,  8p  =  2  x  103  Pa  transports  oxygen  to  the 
membrane  most  effectively  when  the  ratio  q  ~  5-10. 

When  pressure  and  its  drop  are  sufficiently  high,  i.e.,  p2  — 
p2  >  2nq(a/P),  q  =  d/h,  D’Arcy  flow  dominates  over 
diffusion  in  y  direction.  Asymptote  of  catalyst  oxygen  flux 
Jo2  at  high  pressures  and  drops,  i.e.,  p\  —  p2  q2(a/ ft), 
is  independent  of  pressure,  its  drop  and  permeability,  as  the 
diffuser  is  saturated  with  oxygen.  At  lower  pressures,  i.e., 
27tq(a/ p)  <  p\  —  p\  <  q2{a/ P)  the  flux  Jo2  is  proportional 
to  p2  —  p2,  but  it  is  independent  of  the  diffusion  coefficient. 
Note  that  the  flux  appears  to  depend  on  the  dimensionless 
ratio  q  but  not  on  the  geometric  sizes  d  and  h.  When  p2  — 
p2  <  2jtq(a/P),  oxygen  is  transported  by  diffusion  alone. 


4.  Application  of  the  3D  model  to  the  electrode  of  the 
PEFC  with  serpentine  gas  channels 

In  reality,  the  length  of  a  gas  supply  channel  L  in  PEFC 
is  greater  than  the  distance  between  them  d  and  diffuser 
thickness  h:L  d,  L  h.  Then  the  following  relations 
dp/dy  dp/dz ,  V/?||y  hold  everywhere  in  interdigitated 
channel  net,  but  for  the  vicinity  of  the  turn  points  in  con¬ 
tinuous  channel  net.  Contribution  of  the  last  is  small,  re¬ 
spectively,  to  the  rest  of  the  net,  particularly  because  of  the 
low  pressure  drop  between  neighbouring  channels  at  the 
turn  points.  Consequently,  gas  fluxes,  directed  along  the  sup¬ 
ply  channels  in  the  diffuser  are  negligible,  respectively,  to 
those  across  the  diffuser.  Considering  cathode,  the  integra¬ 
tion  along  the  channel  of  the  oxygen  flux  directed  to  the 
catalyst  surface  at  x  =  h  over  the  neighbouring  channels 
results  in  the  following  expression: 


where  E\ ,  E2,  G,  H ,  F  are  given  by  expressions  (7),  where 
pi(z ),  cq(z)  are  the  inlet  channel  pressure  and  oxygen 
concentration,  p2(z),  c\  (z)  are  the  outlet  channel  pressure 
and  oxygen  concentration,  respectively,  a  and  p  are  con¬ 
stants,  given  by  formulae  (6),  d  is  the  distance  between  the 
channels. 

Expression  (9),  applied  to  the  plate  of  serpentine  chan¬ 
nels,  yields  the  total  flux  through  the  membrane  plane  in  the 
following  form: 


Eu  = 
E2i  = 


CQj  -  c\j  exp  {-((Gj  +  Hj)/2)d} 

1  —  exp  {—Gid} 

-CQj  exp  {  —  Gjd}  +  c\j  exp  {-((G/  +  Hj)/2)d} 

1  —  exp  {—Gid} 


(11) 


7 r 

2 h' 


where  E\,  E2 ,  G,  H,  E  are  given  by  formulae  (7),  a , 
P  are  given  by  formulae  (6 ),  d  =  W/(n  —  1),  n  is  the 
number  of  channels,  W  is  the  width  of  the  channel  plate, 
L  the  length  of  one  section  of  the  serpentine  channel, 
pu,  p2i ,  C0i  c\i  are  average  values  of  proper  variables 
in  every  channel,  and  other  notations  are  the  same  as  in 
formulae  (7). 

In  view  of  rather  small  pressure  and  concentration  drops, 
linearisation  can  be  performed  as  follows: 


pu  =  Pin  + 


P  out  P  in 


(12) 


COi  —  Coin  +  (  i  — 


1 


5cq, 


C\i  =  CQi  —  <$Q), 


CO  out  CQ  in 


9 


where  pou t  is  the  outlet  gas  pressure,  pm  is  the  inlet  gas 
pressure,  co0ut  is  the  outlet  reactant  concentration,  coin  is 
the  inlet  reactant  concentration,  n  is  the  number  of  supply 
channels,  h  is  the  thickness  of  the  diffusive  layer,  d  is  the 
distance  between  the  neighbouring  channels,  W  is  the  width 
of  the  channel  plate. 


n°’=ad 


E\{z) 


H(z)  -  G(z) 


H(z )  -G(z)  , 

exp - - - y  + 


E2(z ) 


H(z)  +  G(z) 


H(z)  +  G(z )  ,, 
exp - - - v  dz, 


(9) 
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5.  Current  limitations  in  PEFC  with  2,  3,  4,  5,  7  and  10 
gas  supply  channels  in  serpentine 

5. 1 .  Non  flooded  model 

In  this  section,  the  above  described  model  of  the  gas  trans¬ 
port  in  porous  diffuser  of  the  electrode  of  PEFC  is  applied 
for  the  comparative  analysis  of  the  maximum  current  limita¬ 
tions  of  the  PEFC  for  several  specific  configurations  of  gas 
supply  channels.  The  comparison  is  hold  between  nets  of  2, 
3,  4,  5,  7  and  10  parallel  channels  machined  in  serpentine  on 
a  square  plate.  The  aim  of  the  studies  is  the  identification  of 
most  efficient  in  the  terms  of  the  limiting  current,  channel 
configuration. 

Reduce  the  variety  of  possible  serpentine  configurations  to 
that,  schematically  given  in  Fig.  4,  with  the  effective  average 
pressures  in  channels  derived  from  the  approximation  of 
linear  pressure  drop. 


The  limiting  current  of  the  PEFC  is  determined  by  the 
maximum  oxygen  flow,  converted  on  the  catalyst.  Suppose 
that  the  electrochemical  processes  are  equal  for  both  con¬ 
figurations,  additional  losses  in  the  catalyst  layer  and  those, 
coming  due  to  ohmic  resistance  are  equal  for  all  channel 
configurations.  We  also  neglect  the  effect  of  flooding  of  the 
diffuser  by  liquid  water  on  this  stage.  Set  equal  pressure 
regime  for  both  configurations  of  channels.  The  developed 
model,  applied  for  the  studies  of  air  flows  through  the  dif¬ 
fuser  of  PEFC  over  two,  three  and  10  serpentine  gas  supply 
channels  under  equal  pressures  predicts  higher  values  of  the 
maximum  O2  catalyst  flux  for  smaller  number  of  channels 
in  section. 

Figs.  5  and  6  demonstrate  up  to  15-20%  increase  of  the 
maximum  oxygen  flow  to  the  catalyst,  for  three  and,  particu¬ 
larly,  for  two  channel  serpentine,  respectively,  to  10  channel 
net.  Consequently,  we  expect  the  same  value  for  the  relative 
increase  of  the  limiting  current  of  the  cell. 


inlet 


SERPENTINE  design;  MxN  channels. 

Channel  structure  on  the  plate:  N  sections  x  M  channels. 
25x2, 17x3, 13x4, 10x5,  7x7,  5x10  channels. 

Geometric  parameters  of  the  modelled  device: 

width  of  each  channel  r=  approximately  1mm, 
distance  between  the  channels  d=  approximately  1mm, 
channel  section  length  L=100mm, 
total  plate  width  100mm, 
effective  diffussion  length  h=0.25mm. 


Fig.  4.  Schematic  illustration  of  serpentine  configuration  of  n  parallel  gas  channels  in  a  PEFC. 


EelAt-ive  oxygen  flow  increment-  of  t  llchAnnel.'  a qii rut- 

10  ||  charnel.;  iri  SERPENTINE  5ect-ioru  on  10X10  cm  plAt-e 
for  irilet-  Air  prej-jure  Arid  pre^jijre  drop 


Fig.  5.  Relative  increment  of  the  maximum  oxygen  flow  in  diffuser  over  two  channel  serpentine,  respectively,  to  10  channel  serpentine  for  inlet  pressure 
pi  and  pressure  drop  A p. 
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Eelat-ive  oxygen  flow  inc remerit-  of  3  llchumelj  igiiri5t - 

10  ||  channel  5  in  SERPENTINE  ject  iorij  on  10X10  cm  pi  it? 

for  inlet-  iir  pressure  and  pressure  drop 


Fig.  6.  Relative  increment  of  the  maximum  oxygen  flow  in  diffuser  over  three  channel  serpentine,  respectively,  to  10  channel  serpentine  for  inlet  pressure 
pi  and  pressure  drop  A p. 


However,  more  realistic  situation  can  be  modelled  in  the 
following  way.  Fix  the  inlet  pressure,  and  vary  the  oxygen 
flow  in  the  gas  supply  channel,  following  the  electric  cur¬ 
rent  of  the  cell.  Then,  A p  will  be  specific  for  each  channel 
configuration,  i.e.  two,  three,  etc.  channels  in  a  serpentine 
section,  for  the  same  current.  This  predictions  of  the  model, 
based  on  the  evaluation  for  the  airflow  values  approximately 
equal  to  40  cm3  min-1  A-1  result  in  13%  limiting  oxygen 
flux  increase  for  three  channel  serpentine,  respectively  to  10 
channel  design. 

The  predicted  values  of  the  relative  increment  of  the  lim¬ 
iting  current  (due  to  oxygen  transport)  for  various  numbers 
of  channels  in  serpentine  are  given  in  Table  1. 

5.2.  Flooded  model 


and  obtain  precise  information  about  the  quantities  of  liquid 
water  in  various  points  of  the  diffuser,  we  suppose  that  the 
water  reduces  the  volume  of  the  porous  in  the  diffuser  by 
maximum  a  quarter,  whereas  the  removal  of  the  water  by 
convection  is  complete. 

Consider  the  following  flooded  model,  where  12%  of 
porous  in  the  diffuser  over  the  channels  and  25%  or  porous 
in  the  diffuser  between  the  channels  of  each  serpentine  sec¬ 
tion  are  flooded  by  liquid  water.  Then  the  porosity  in  each 
domain  is  equal  to  £\  =  0.8,  £2  =  0.7,  £3  =  0.6,  respec¬ 
tively.  Account  for  the  effective  tortuosity  and  porosity  of 
the  diffuser  will  be  carried  out  through  the  Bruggeman-type 
correction,  which  estimates  effective  oxygen  diffusion  coef¬ 
ficient  with  the  following  relation: 


OZ  =  D0,  si 


O2I 


2^1 


(13) 


In  order  to  take  into  account  the  presence  of  liquid  wa¬ 
ter  in  the  porous  of  the  gas  diffuser,  the  effective  porosity 
will  be  introduced  now.  It  describes  the  porous  media  when 
the  porous  are  partially  filled  with  liquid  water  in  thermody¬ 
namic  and  hydrodynamic  equilibrium.  The  quantity  of  liq¬ 
uid  water,  retained  by  gas  diffuser,  depends  on  the  nature 
of  the  mass  transport  in  it.  Due  to  the  pressure  gradient  be¬ 
tween  the  neighbouring  serpentine  sections,  forced  convec¬ 
tion  helps  removing  the  majority  of  water  from  the  diffuser 
there.  Where  the  gas  channels  lay,  a  certain  part  of  the  liq¬ 
uid  is  removed  from  the  diffuser  just  over  them  due  to  the 
channel  gas  flow.  The  diffuser  over  the  shoulders  between 
the  channels  within  a  section  of  the  serpentine  is  probable 
to  be  mostly  flooded,  since  the  water  does  not  have  an  im¬ 
mediate  exit  to  the  gas  channel  and  neither  the  convection 
contributes  to  its  removal  due  to  the  absence  of  the  pressure 
drop  between  such  channels.  Although,  it  is  hard  to  measure 


where  Do2  is  determined  by  formula  (3),  £/  and  r  are  poros¬ 
ity  and  tortuosity  factors  of  the  diffuser  in  i  domain.  Three 
different  domains  of  the  diffuser  are  those  over  the  channel, 
between  the  channels  within  a  section  and  between  the  sec¬ 
tions  of  the  serpentine.  The  above  discussed  modelling  of 
the  diffuser  with  the  number  of  channels  in  serpentine  vary¬ 
ing  from  two  to  ten  results  in  the  following  predictions  of 
the  oxygen  transport  limitations  of  the  cell  for  two  pressure 
regimes,  given  in  Tables  2  and  3. 

Comparing  Tables  1  and  2  we  note  that  with  the  account 
for  the  liquid  water,  present  in  the  porous,  the  efficiency 
of  the  serpentine  with  low  number  of  channels  in  sections 
becomes  even  more  evident.  The  limiting  oxygen  flow  in¬ 
crease  for  M  =  3  channels  in  serpentine  now  equals  to 
22-25%  against  13%  predicted  with  no  account  for  water 
flooding.  The  pressure  regime  is  chosen  so  that  the  airflow 
in  the  channels  is  independent  on  the  channel  configuration; 
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Table  1 


Relative  increment  of  the  limiting  current,  coming  due  to  oxygen  transport  limitations  in  N  channel  serpentine,  respectively,  to  10  channel  serpentine  for 
equal  air  flow  to  electric  current  ratio 


M  channels  x  N  sections 

10  x  5 

7x7 

5  x  10 

4  x  13 

3  x  17 

2  x  25 

(Pi  —  1.6  atm) 

(Pi  =  1.6  atm) 

(Pi  =  1.6  atm) 

(Pi  =1.7  atm) 

(Pi  =1.7  atm) 

(Pi  =  1.75  atm) 

Maximum  current  /Mch-floch  (%) 

0 

3 

6 

9 

13 

21 

Pressure  drop  A  P  =  P-m  —  Pout  (atm) 

0.05 

0.1 

0.2 

0.33 

0.57 

1.25 

Table  2 


Relative  increment  of  the  limiting  current,  coming  due  to  oxygen  transport  limitations  in  partially  flooded  model  in  N  channel  serpentine,  respectively, 
to  10  channel  serpentine  for  constant  air  flow  to  electric  current  ratio  for  A P  =  0.075  atm  for  10  channel  geometry 


M  channels  x  N  sections 

10  x  5 

7x7 

5  x  10 

4  x  13 

3  x  17 

(Pi  -  1.6  atm) 

(Pi  —  1.6  atm) 

(Pi  -  1.6  atm) 

(Pi  =1.7  atm) 

(Pi  =1.7  atm) 

Maximum  current 

0 

4 

10 

15 

24 

4lch— flOch  (%) 

Pressure  drop 

0.075 

0.14 

0.27 

0.46 

0.85 

AP  -  Pin  -  Pout  (atm) 

Table  3 

Relative  increment  of  the  limiting  current,  coming  due  to  oxygen  transport  limitations  in  partially  flooded  model  in  N  channel  serpentine,  respectively, 

to  10  channel  serpentine  for  constant  air  flow  to  electric  current  ratio  for  A  P  = 

=  0.035  atm  for  10  channel  geometry 

M  channels  x  N  sections 

10  x  5 

7x7  5  x 

10  4  x  13 

3  x  17 

2  x  25 

(Pi  —  1.6  atm) 

(Pi  =  1.6  atm)  (Pi  : 

=  1.6  atm)  (Pi  —  1.7  atm) 

(Pi  =1.7  atm) 

(Pi  -  1.7  atm) 

Maximum  current 

0 

4  9 

15 

23 

39 

4/Ich-fl0ch  (%) 

Pressure  drop 

0.035 

0.06  0.13 

0.21 

0.4 

0.9 

AP  =  Pin  Pout  (atm) 

Table  4 

Relative  increment  of  the  limiting  current,  coming  due  to  oxygen  transport  limitations  in  dynamical  flooded  model  in  N  channel  serpentine,  respectively, 

to  10  channel  serpentine  for  constant  air  flow  to  electric  current  ratio  for  A  P  = 

=  0.075  atm  for  10  channel  geometry 

M  channels  x  N  sections 

10  x  5 

7x7 

5  x  10 

4  x  13 

3  x  17 

(Pi  —  1.6  atm) 

(Pi  —  1.6  atm) 

(Pi  -  1.6  atm) 

(Pi  =1.7  atm) 

(Pi  =1.7  atm) 

Maximum  current 

0 

4 

10 

16 

25 

4lch-fl0ch  (%) 

Pressure  drop 

0.075 

0.14 

0.27 

0.46 

0.85 

AP  =  Pin  -  Pout  (atm) 

for  example,  the  values  in  Tables  3  and  5  are  calculated 
for  the  airflow  ^50  cm3  min-1  A-1.  Variations  in  the  air 
flow  and  consequently  in  the  pressure  regime,  reported 
in  Tables  2  and  3  effects  the  predictions  for  the  maxi¬ 
mum  catalyst  oxygen  flux  only  slightly.  However,  it  may 
be  difficult  to  maintain  sufficient  air  supply,  independent 
on  different  channel  sets,  because  the  required  pressure 
drop  should  be  high  for  low  number  of  channels  in  a 
section  M. 


Besides  the  above  described  modelling  of  the  electrode 
flooding,  the  dynamic  model  of  the  flooding  effect  in  the 
electrode  is  developed.  It  is  presumed  that  the  percentage 
of  the  porous,  available  for  gas  transport  depends  on  the 
pressure  gradient  S/p  exponentially,  so  the  effective  porosity 
varies  from  £\  =  0.8  when  all  the  water  is  removed  by 
the  shear  flow  with  V p  =  0.2  atm  cm-1  pressure  gradient, 
to  s\  =  0.55  when  S/ p  =  0.  Such  approach  results  in  the 
following  predictions,  given  in  Tables  4  and  5: 


Table  5 


Relative  increment  of  the  limiting  current,  coming  due  to  oxygen  transport  limitations  in  dynamical  flooded  model  in  N  channel  serpentine,  respectively, 
to  10  channel  serpentine  for  constant  air  flow  to  electric  current  ratio  for  A P  =  0.035  atm  for  10  channel  geometry 


M  channels  x  N  sections 

10  x  5 

(Pi  —  1.6  atm) 

7x7 

(Pi  —  1.6  atm) 

5  x  10 

(Pi  —  1.6  atm) 

4  x  13 

(Pi  —  1.7  atm) 

3  x  17 

(Pi  —  1.7  atm) 

2  x  25 
(Pi  -  1.7  atm) 

Maximum  current 

4/Ich-fl0ch  (%) 

0 

3 

8 

13 

22 

40 

Pressure  drop 

AP  =  Pin  Pout  (atm) 

0.035 

0.06 

0.13 

0.21 

0.4 

0.9 
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The  predictions  of  dynamical  and  partially  flooded  models 
appear  to  be  similar,  as  the  comparison  of  Table  2  with 
Tables  3-5  evidences.  For  low  pressure  drops,  dynamic  flood 
model  predicts  slightly  lower  performance  of  seven,  five 
and  four  channel  serpentines,  respectively,  to  the  predictions 
of  the  partially  flooded  model — 3,  8  and  13%  against  4,  9 
and  14%  increase  (see  Tables  3  and  5),  compared  with  10 
channel  serpentine.  At  higher  pressure  drops  the  dynamic 
model  shows  slightly  higher  results  -  16%  for  M  =  4  and 
25%  for  M  =  3,  against  15  and  24%,  respectively,  for  four 
and  three  channels  in  partially  flooded  model  (see  Tables  2 
and  4,  also  compare  the  results  for  M  =  2  in  Tables  3  and  5). 

6.  Discussion  and  conclusions 

The  experimental  measurements  demonstrated  ^25% 
higher  cell  current  for  three-channel  parallel  serpentine 
against  10  channel  configuration  in  the  high  current  regime, 
when  the  air  flow  in  air  supply  channels  followed  the  pro¬ 
duced  electric  current  at  52  cm3  min-1  A-1  (the  maximum 
current  was  limited  at  0.6  A  cm-2  due  to  technical  reasons). 
This  experimental  data  are  in  agreement  with  predicted 
theoretical  values  and  demonstrate  the  same  tendency,  i.e. 
the  efficiency  of  the  PEFC  at  high  currents  increases  with 
the  decrease  of  the  number  of  supply  channels  M  in  a  ser¬ 
pentine  section,  provided  the  air  flux  in  supply  channels  is 
maintained  independent  on  M.  However,  it  requires  high 
pressure  drops  A p  for  serpentines  with  M  =  2  and  M  —  3. 
Developed  with  no  account  for  liquid  water  model  predicts 
13%  current  increase  for  three-channel  net,  respectively,  to 
10  channel  net  (see  Table  1).  It  agrees  with  the  experiment 
only  qualitatively.  The  difference  between  experimental  and 
model  values  can  be  explained  by  the  effect  of  liquid  wa¬ 
ter,  present  on  the  cathode  side  in  the  porous  diffuser,  but 
very  difficult  to  model.  Due  to  the  big  number  of  sections 
and  higher  A p  between  the  sections  (which  is  necessary  to 
provide  fixed  air  supply  (cm3  min-1  A-1))  for  M  =  2  and 
M  =  3,  respectively,  to  M  =  10  channel  serpentines,  the 
air  is  forced  to  flow  in  the  diffuser  between  the  channels  of 
neighbouring  sections.  This  effect  helps  removing  the  liquid 
water  from  the  porous,  since  molecular  diffusion  becomes 
combined  with  forced  convection  in  the  matter. 

Really,  developed  with  the  account  for  the  liquid  water, 
models  predict  22-25%  higher  performance  of  three-channel 
serpentine,  respectively,  to  10  channel  serpentine  in  the  terms 
of  the  limiting  current.  These  values  are  in  agreement  with 
experimental  measurements.  The  predictions  of  dynamical 
and  partially  flooded  models  appear  to  be  similar,  as  we 
can  see  that  from  the  comparison  of  Tables  2  and  3  with 
Tables  4  and  5,  respectively.  Dynamic  flood  modelling  yields 
some  weaker  advantage  for  seven,  five  and  four-channel 
serpentines  than  partially  flooded  model  does  at  low  pressure 
drops  (see  Tables  3  and  5). 

However,  the  higher  the  pressure,  the  better  the  liquid  wa¬ 
ter  is  removed  from  the  pores  and  the  higher  the  predictions 


of  the  dynamic  model  become.  It  also  predicts  higher  limit¬ 
ing  currents  than  it  is  expected  with  constant  flooding  (see 
Tables  2  and  4  for  M  =  4  and  M  =  3,  also  see  the  results 
for  M  =  2  in  Tables  3  and  5).  Although  the  difference  be¬ 
tween  the  predictions  in  Tables  4  and  5  is  about  of  2-3%  for 
three,  four  and  five  channels,  this  important  dependence  on 
the  pressure  regime,  accounted  for  in  the  dynamical  model, 
is  in  agreement  with  few  available  experimental  observa¬ 
tions,  which  demonstrate  that  the  limiting  current  generally 
increases  slightly  with  the  rise  of  the  pressure  and  its  drop. 

The  agreement  between  theoretical  predictions  of  the  de¬ 
veloped  model  and  the  data  of  experimental  measurements 
indicates  on  the  correct  account  for  the  complex  of  the  phys¬ 
ical  phenomena,  responsible  for  the  differences  in  current 
limits  between  geometric  configurations  and  its  adequate 
modelling. 

It  is  important  to  underline  that  maximum  current  increase 
for  two  channel  net  hchlhoch  could  reach  35-40%,  but  to 
obtain  equal  air  supply  flux  for  M  =  10  and  M  =  2  chan¬ 
nel  geometry,  the  pressure  drop  in  the  last  case  should  be 
A  P2  channel  net  ^  1.8  atm  if  A  Pi 0  channel  net  ^  0.075  atm  and 
A  P2  channel  net  ^0.85  atm  if  A  Pi 0  channel  net  ^  0.035  atm,  re¬ 
quiring  high  inlet  pressure.  Due  to  technical  reasons,  such 
a  regime  was  impossible  to  set  and  the  inlet  pressure  was 
fixed  at  ^1.6  atm. 

The  results  obtained  and  discussed  above  lead  to  the  fol¬ 
lowing  general  conclusions: 

•  Due  to  the  forced  convection  mechanism,  serpentine  gas 
supply  channel  nets  with  small  number  of  channels  M  in 
each  section  provide  higher  limits  on  directed  to  catalyst 
O2  flux  in  diffuser  and  therefore  higher  limits  on  the  max¬ 
imum  current.  Gas  supply  configurations  with  two  and 
three  channels  were  identified  as  the  best  for  the  high  lim¬ 
iting  current.  (Dependently  on  the  pressure  regime,  up  to 
25%  increase  for  3  channel  net  and  up  to  35-40%  in¬ 
crease  for  two  channel,  respectively,  to  the  serpentine  of 
M  =  10  channels  in  a  section) 

•  Negative  effect  of  water  flooding  in  diffuser  is  reduced 
when  serpentines  with  low  number  of  channels  in  a  section 
M  and  high  number  of  sections  N  are  employed. 

•  Very  small  number  of  channels  (M  =  1,  2)  in  a  serpentine 
section  may  result  in  significant  deplete  of  oxygen  along 
the  channels  due  to  insufficient  air  supply  at  low  pressure 
drops  and  require  high  pressure  drops. 

Modelling  oxygen  transport  in  the  diffuser  of  the  PEFC 
with  pressure  1  atm  <  P  <  2  atm,  Oatm  <  A P  <  0.5  atm 
and  the  number  of  channels  in  each  serpentine  section  M  = 
2,  3, . . .  ,  10,  we  conclude  that  serpentine  channel  nets  with 
M  =  3  and  M  =  2  channels  in  each  section  provide  higher 
limits  on  catalyst  oxygen  flux  in  diffuser  and  therefore  es¬ 
tablish  higher  limits  on  the  maximum  current:  up  to  13% 
more  for  N  =  3  and  20%  increase  for  N  =  2,  respectively, 
to  M  =  10  channels  in  serpentine.  (See  Figs.  5  and  6). 

Modelling  oxygen  transport  in  the  diffuser  of  the  PEFC 
with  liquid  water  flooding  at  the  air  flow  ^50  sec  min-1  A-1 , 
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and  the  number  of  channels  in  a  section  M  =  2,  3,  . . .  ,10, 
we  found  the  predictions  of  the  model  in  agreement  with 
experimental  data. 

•  With  the  account  for  flooding  of  the  diffuser  by  liquid  wa¬ 
ter,  serpentines  with  M  =  4,  3  or  2  channels  in  a  section 
appear  to  be  more  effective  than  10  channel  design,  pro¬ 
viding  even  higher  limits  on  the  maximum  current  that 
the  model  without  flooding:  up  to  25%  for  M  =  3  (see 
Tables  2  and  4).  Serpentine  with  M  =  2  would  provide 
35-40%  higher  current  limits  than  M  =  10  channel  net 
(see  Tables  3  and  5),  but  it  would  also  require  high  pres¬ 
sure  drops  to  maintain  required  air  flux. 

•  Dynamic  modelling  of  flooding,  accounting  for  liquid  wa¬ 
ter  removal  by  shear  gas  flow  in  diffuser,  dependent  on  the 
pressure  gradient,  predicts  few  percent  better  performance 
of  the  PEFC  working  with  higher  pressure  and  pressure 
drop.  This  agrees  with  experimental  observations. 

The  effect  of  flooding  of  the  electrode  can  be  investigated 
deeper  in  future  studies  together  with  further  adjustments  of 
the  model.  Additional  experimental  tests  are  foreseen. 
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